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Abstract 

We unify all existing results on the change of the speed of low-energy photons due 
to modifications of the vacuum, finding that it is given by a universal constant times the 
quotient of the difference of energy densities between the usual and modified vacua over 
the mass of the electron to the fourth power. Whether photons move faster or slower than 
c depends only on the lower or higher energy density of the modified vacuum, respectively. 
Physically, a higher energy density is characterized by the presence of additional particles 
(real or virtual) in the vacuum whereas a lower one stems from the absence of some virtual 
modes. We then carry out a systematic study of the speed of propagation of massless 
particles for several field theories up to two loops on a thermal vacuum. Only low-energy 
massless particles corresponding to a massive theory show genuine modifications of their 
speed while remaining massless. All other modifications are mass-related, or running mass- 
related. We also develop a formalism for the Casimir vacuum which parallels the thermal 
one and check that photons travel faster than c between plates. 
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1. Review of existing results 



The usual vacuum in a Quantum Field Theory is a state characterized by the absence 
of real particles and classical fields and by its Minkowskian geometry. Electromagnetic and 
gravitational fields as well as massless particles propagate though it with the same constant, 
Lorentz invariant speed, c. When the vacuum is modified, so is the speed of propagation of 
particles and fields. This quantum field theoretical effect has been analyzed within QED 
for low-energy photons in several cases that we now review. 

Adler first studied the propagation of photons in an anisotropic vacuum given by an 
external constant uniform magnetic field, B, and obtained the result (h = c = eo = 1) [1] 

8 r, .2/1 

vu = l-—a — jsm d <1, 
45 mj 

v± = 1 — J sm"' 6 < 1 , 

45 

for polarizations respectively coplanar with and perpendicular to the plane defined by 
B and the direction of propagation. Here, me is the mass of the electron, a is the fine 
structure constant and 6 is the angle between B and the direction of propagation. 

For a (homogeneous and isotropic) Robertson- Walker gravitational background with 
Friedmann cosmology, Drummond and Hathrell obtained [2] 



v = l + ^aGN^ >1, (1.2) 
45 mi 



where Gn is Newton's constant and p and p are the energy density and pressure (which 
satisfy by the energy condition [3], p+p > 0). Notice that the speed of light is now greater 
than c. 

For a homogeneous and isotropic thermal vacuum, one of us first computed the change 
in the speed of light [4] . That result was corrected in Ref. [5] and is reobtained here in 
Sect. 5 with the result 

v = l a^—r <1, 1.3) 

2025 mi ' ^ ^ 

where T is the temperature and kB = 1- 

Finally, Scharnhorst [6] and, independently. Barton [5] worked out the propagation of 
light in a (anisotropic and non-translational invariant) Casimir vacuum with the result 

V ^ 1 + >1, (1.4) 

8100 a'^mj ' ^ ' 

where a is the distance between plates and 6 is the angle between the direction of prop- 
agation and the normal to the plates. Again, photons propagate faster than c in this 
case. 
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There are several features worth mentioning in these four expressions. First, we note 
that the electron mass plays the role of an infrared regulator, so that the result is analytic 
for low energy and momenta. Also, mg provides the natural dimensions to compensate 
for those of the modified vacuum and, thus, sets the scale for this effect to be physically 
observable. All the above results are perturbative and should only be believed when small. 

— * 

The problem is that they are extremely small since any easily accessible values of B, 
Gn{p -\- p)i T or are negligible when compared to rUe. In other words, electrons are 
too heavy in this context. 

A second relevant observation is that the modification of the polarization-averaged 
speed of light is always proportional to the numerical factor 11. This was already pointed 
out with respect to Eqs. (1.3), (1.4) in Ref. [5], where the author noted that this factor is 
common to all coherent light-by-light scattering [7]. Notice, however, that the same factor 
appears in Eq. (1.2) which is unrelated to light-by-light scattering, hinting at some deeper 
interpretation of the above results. 

It is, nevertheless, hardly arguable that the most enticing feature of the above results 
is the sign of the modification of the speed of light: some vacua lead to ^; < c, some to 
V > c. This last possibility has motivated further study [5] [8] [9] [10]. The reader will find 
there a controversial discussion of phase, group and causal signals and, in particular, of 
the effects of > c on causality. Let us just say that, in accordance with ref. [2], we find 
no grounds for violation of microcausality. 

The purpose of our work is to deepen and complete our understanding of the mod- 
ifications of the speed of light. We do this by proceeding from the more physical to the 
more field-theoretical aspects of the problem. We first study the already known results 
presented above. We will give a unified and very physical interpretation of them which 
explains in an extremely simple manner the sign of the variation of the speed of light. 
We then concentrate on the thermal vacuum and consider the modification of the speed 
of (high and low energy) massless particles for different theories up to two loops. This 
will require a precise definition of what is meant by speed of massless particles which is 
presented in Sect. 3. A comprehensive one-loop study shows that field theories fall into 
two categories for the purpose of this study: those with an explicit mass, m, and those 
with an intrinsic mass scale, A. This is done in Sect. 4. In Sect. 5, the two-loop study 
for a theory of the first category, QED, is performed. This reproduces Eq. (1.3) at low 
energies and leads to a new result at high energies. The two-loop study for a theory of the 
second category, massless gcf)'^, is performed in Sect. 6. This is mainly of field theoretical 
interest as in the Standard Model only the quarkless QCD sector has no explicit mass. 
Confinement only allows massless particles beyond the deconfinement temperature which 
does not correspond to the low temperature cold heat bath situation we are considering 
here (T <^ m or T <^ A). Still, it is not excluded that these results will be of some 
relevance to photon propagation if almost massless pseudogoldstone bosons, like invisible 
axions, exist. In Sect. 7, a comparison of the thermal and Casimir vacua is performed. 
This allows us, first, to better understand the meaning of a modified vacuum and, second. 
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to introduce a unified formafism which explains, both mathematically and theoretically, 
the different signs of the modifications of the speed of light. We finally summarize our 
most relevant results. 



2. The unified formula 



We mentioned in the previous section that all the polarization-averaged modifications 
of the speed of light are proportional to the numerical factor 11 and inversely proportional 
to mg. This hints at the possibility of writing them in a unified form, which thus digs out 
the common phenomenon underlying the physics behind them. This is indeed so, as we 
now show. 

Consider first the magnetic field vacuum. From Eq. (1.1) one finds the averaged speed 
of light 

^=\£d0 sin^ (^11 (^) + v^i9)) = 1 - ^ ^ . (2-1) 
which, recalling that the energy density for a magnetic field is given by 

PB = ^>0, (2.2) 



can be written as 

44 2 PB 

a — 

135 m- 



v = l-^ a"" ^ . (2.3) 



For the Friedmann-Robertson- Walker vacuum, we recall that a radiation dominated 
universe, as corresponds to T <C me, is described by p — |. Moreover, the gravitational 
energy density is negative 

PG = -p<0, (2.4) 



so that Eq. (1.2) can be written as 



^ = 1 - ^ « (GNml) ^ . (2.5) 



e 



On its own, the cold temperature vacuum is characterized by its Planck energy density 



so that Eq. (1.3) reads 



PT = ^T'>0, (2.6) 



" = '-135" ^^-'^ 
4 



Finally, the averaged speed of light for the Casimir vacuum is, from Eq. (1.4), 

v=- [ de sm9v(9) = l + ^^ ( ■ (2.8) 
2 Jq ^ ^ 24300 \amej ^ ^ 

Recalling the energy density of the Casimir vacuum [11], 
the speed of light can be written as 

^; = 1 _ ii ^2 . (2.10) 

135 mi ^ ' 

Note that Eqs. (2.3), (2.7) and (2.10) are identical, and (2.5) also coincides with them 
if we substitute GNm\ by a. We therefore find the following remarkable result: the low 
energy modification of the speed of light is proportional to the ratio of the energy density 
of the modified vacuum (as compared to the usual one) over mg, with a universal numerical 
coefficient and the corresponding coupling constants. If the modified vacuum has a higher 
energy density, then v < c. If it is lower, then v > c. This is what explains the sign of the 
change of the speed of light. 

The general equation we have found for the speed of light in a modified vacuum allows 
us to make further predictions. As an instance, low-energy photons should travel in an 
external electric field, E, with a speed 

y = l-^a^^. (2.11) 
135 mi ^ ' 

Another remarkable consequence of the general equation we have found is the fact that only 
for an asymptotically flat universe would the speed of light remain c. Matter, radiation 
and gravitational energy add to zero, which is the relevant physical quantity that affects 
the velocity of low-energy photons at this order of perturbation theory. 

Let us finish this section mentioning one surprising feature. Notice that the numerical 
coefficient is the same for the gravitational vacuum as for the others, in spite of its being 
unrelated to light-by-light scattering. Furthermore this factor 11 also appears in the co- 
efficient of the Euler-Poincare characteristic spin i contribution to the gravitational trace 
anomaly [11] [12]. These connections will be discussed elsewhere. 



3. Definitions and framework 

Let us first define what we mean by speed of light, or better, speed of a massless parti- 
cle. Consider a massless scalar particle which moves in a thermal vacuum. Its propagator 
is given by 

^ n (0) = , (3.1) 
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where Z is the field renormahzation constant, 11 the zero temperature self-energy and its 
temperature correction. The separation of and in 11^ reflects the Lorentz-symmetry 
breaking induced by the heat bath. 

We will define the speed of a massless particle for real, rather than virtual, particles. 
This is because for gauge theories definitions on the pole, even at finite temperature, ensure 
gauge invariance [13]. The energy-momentum dispersion relation which characterizes the 
propagation of a real particle is obtained from imposing that the denominator of Eq. (3.1) 
vanishes. This leads to the implicit solution go = ^{fl)-, 9 = l^l- At one-loop order one 
finds an explicit solution independent of 11 , 

E2(,)=g2 + ^«(,^g2). (3.2) 



Note that the standard vacuum polarization, H, does not break Lorentz invariance and, 
thus, does not change the energy-momentum dispersion relation. 

In general, perturbation theory allows to obtain an explicit solution 

E^{q) = q^^j{q,T,m) , (3.3) 



where /(g, T, m) depends on 11 and 11^. Its dependence on momentum, temperature and 
virtual particle mass has been specified. If the theory is massless, m = 0, an intrinsic mass 
scale. A, will appear beyond one loop through the running of the coupling constant. Then, 
/(g, T, m) should be understood as /(qf, T, A). Certainly the dependencies on m and A are 
different. 

The speed of the massless particle is then given (as for photons traveling through a 
medium that changes the relation between E and p) by 

Expansions of potential relevance are 



(/(?, T, m) - g/'(g, T, m)) , » /, g » /' , (3.5) 



with ~ 1, 



2v/(g,r, m) 



and 

V , /'(g, T, m) « g, / » g2 , (3.7) 

with V ^ 0, in spite of going, to one loop, as the inverse of the square root of the coupling 
constant. 
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Three particular cases will be of interest. First, let / be constant in 



f{q,T,m) = X\T,m) . (3.8) 

The whole effect of the heat bath is, then, to give a mass to the particle, which however 
remains free, see Eq. (2.3). The speed of the particle is modified as corresponds to 
acquiring a mass: from Eq. (3.5) and Eq. (3.7), as Eq. (3.6) never applies. 



vc^l , 9>A (a) 



^ , A»g (b) 



(3.9) 



A(T, m) 

Second, let us consider the case when no mass is induced at low energies, 

lim/(g,r,m) = 0. (3.10) 

q-»0 

Then, v{q) is a genuine modification of the speed of the massless particle at low energies, 
since it remains massless. The dispersion relation is usually non-free. If, furthermore, 

/(g,T,m)-0(g2) , g«r (3.11) 
(recall that always T <^ m), then Eq. (3.5) always holds in perturbation theory and 

v{q<^T)c^l + i^^^^^^^ , r«l. (3.12) 

This is the situation which corresponds to the results reviewed in the introduction. In the 
particular case f{q,T,m) ~ q^ for all momenta the dispersion relation would be free but 
the speed of the massless particle would have changed according to Eq. (3.12) for all values 
of q, i.e. it would be a different constant now. We have not encountered this situation in 
any of the theories we have investigated. 

Third, let us consider the case when no mass is induced at high energies, 

lim /(?, T, m) = . (3.13) 

g— >oo 

Then, at high energies Eq. (3.5) holds and represents a genuine modification of the speed 
of the massless particle at high energies, for it remains massless. The dispersion relation 
is never free. 

We have been implicitly assuming that f{q,T,m) is real. This is, however, usually 
not so [14], but imaginary contributions of the self energy are not relevant to the issue of 
the speed of the massless particle. Also, no zero momentum ambiguities [15] [16] appear in 
what follows. 
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For massless Dirac fermions the propagator is given by 

go7o -q-1- S(Qo - -Q-l)- Ss,t(?0' 9^)?o7o + S„,t(Qo' 7 ' ^^'^"^^ 

with a notation similar to the used in Eq. (3.1). Notice that no temperature contribution 
of other Dirac matrix structures have been included, because they will not appear for the 
theories and perturbative order considered in this work. The dispersion relation is given 

by 

q', (1 - S(g2 - g2) _ S,,^(go^g2))2 ^ ^2 (1 _ ^.^^2 _ ^2) _ S,,^ (^^^ g^)) 2 _ ^^^^^^ 

At one-loop order this leads to 

E\q) = q'{l + 2 S« (g^ q^) - 2 s(^i(g^ q^)) (3.16) 
or, with the notation of Eq. (3.3), to 

/(?, T, m) = 2q' (s^ q') - ^i^k?', q')) ■ (3-17) 

For massless gauge bosons the corresponding analysis is slightly more complicated [4] . 
The result to one-loop order is 

E\q) = q'- lis'' - m ^«^(^^ q') , (3.18) 

where is the spatial unit vector and Ilij,T are the spatial components of the temperature 
contribution to the self-energy tensor n^^, so that 

f{q,T,m) = -^iS^^-q^q^) U^^^^{q',q^) . (3.19) 

The one-loop expressions also hold at, say, two loops when there are no one-loop 
temperature corrections. This will happen for our two-loop QED study, because we will 
always neglect the thermal contribution due to massive virtual particles which are always 
exponentially suppressed by a huge Boltzmann factor. This is because we consider, as said 
in the introduction, only cold heat baths, T <^ m. 

Summing up: our test particle is real and massless and we only thermalize massless 
virtual particles, so that if the theory has a massive particle it only enters in the form of 
a virtual non-thermalized contribution. 



4. One— loop results 

Let us consider the propagation of massless particles in a cold heat bath, using pertur- 
bation theory at one loop. We are here not concerned with hot heat baths, which require 
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a different formalism [17], because our initial motivation was the smallness of the effects 
reviewed in the introduction, and their extension to high energies and massless theories, 
but without changing the conditions defining the vacuum. We will use the real time for- 
malism [18], which at the one-loop level reduces to the original one-component formalism 
[19]. Thus, a massless scalar boson propagator reads 



DT{k,k') = 5^^\k-k') 



r/r9\ 

+ ^;^ ) 



e = kl- P (4.1) 



+ irj _ I 

while a gauge boson propagator in the Feynman gauge is given by 

D>^''{k,k') = -g'^-DT{k,k') (4.2) 
For a massless Dirac fermion one has 

i 2n 



5{k^) 



(4.3) 



with ^ = kij,^^. The above expressions separate the zero temperature from the finite 
temperature contributions. A heat bath has, of course, other effects on the propagation 
of a particle (e.g. diffusion), but the ones which cannot be told from the virtual particle 
effects, because they conserve energy and momentum, are included in the modified vacuum 
and arc given by Eq. (4.1,2,3). 

Consider first massless (me = 0) QED. For the photon one finds from Eq. (3.19) and 
(4.3) that f{q,T) is constant, i.e. Eq. (3.8) applies and 

A^(T) = 2-^ , (4.4) 

where we have stressed that any of the two virtual electrons thermalizes, and e is the 
electron charge. For the massless electron, from Eq. (3.17), (4.2) and (4.3), again f{q,T) 
comes out constant 

where we have separated the thermalized photon from the thermalized electron contribu- 
tions. The constancy of f{q,T) at one loop turns out to be a general one-loop result for 
theories without explicit masses. Indeed, in SU (Nc) Yang-Mills theory without quarks one 
obtains (ghosts are thermalized as scalar bosons) 

Aj(r) = nJ-^ , (4.6) 

where Nc is the number of colors and gs the strong coupling constant. Also for massless 
|y^^, where the tadpole is the only diagram which appears at one loop, the result is 
constant 

^UT) = ^ ■ (4.7) 
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Consider now the other type of theories we are interested in: those with a virtual 
non-thermahzed massive particle. At one-loop order there are no thermal contributions in 
QED, as the test particle, being massless, has to be the photon and the virtual ones are all 
massive electrons. This will be important when QED is studied at two loops in the next 
section. The Weinberg-Salam model has, however, diagrams which contribute, specifically 
the one with neutrinos both as test and thermalized particles and the Z° as the virtual 
massive particle. The result is 



^ + -r4- / dx -In 



6 47r2g Jq + 1 



M| - 4qTx 



M2 + AqTx 



(4.8) 



where ow = „ . ^ ^ 3—, 6w is Weinberg's angle and is the Z° mass. This satisfies 

Eq. (3.10) but not Eq. (3.13). Thus, the neutrino remains massless at low energies and 
its low energy velocity is modified according to Eq. (3.12), 

v{qT«M!) = l-g^^'4(^y + o('^] , (4.9) 



45 \MJ \ M| 



which is always a small correction for cold heat baths, while at high energies it goes like 



1 fT\^ /M,4 



^;(9T»M|) = 1-^^^(^^J +O^^J , T<g, (4.10) 

where we have specified the condition for perturbation theory to hold. Notice that Eq. 
(4.10) is like Eq. (3.9a) so that the neutrino has indeed acquired a mass at high energies. 
This can also be seen from the fact that Eq. (4.10) allows the limit Mz — to be taken, 
without being modified. But, as analyzed previously in this section, for — the whole 
temperature effect should give the neutrino a mass. Indeed, putting = in Eq. (4.8) 
gives (recall that still only neutrinos are thermalized) 

A^(r)=^^^, (4.11) 

so that Eq. (4.10) can be written as 

^;(gT » Ml) = 1 - ^Ip^ , A.«g, (4.12) 

exactly as Eq. (3.9a). Thus, the neutrino does not acquire a mass at low energies, but it 
does so at high energies, and the mass it acquires is precisely the one it would have gotten 
in the massless theory. In other words 

hm f{q, T, M,) = f{q, T, 0) = (T) (4.13) 

q— >oo 
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The result Eq. (4.9) (with v < 1 and its quartic T dependence), the result Eq. (4.10) with 
V < 1 and its quadratic T dependence) and its interpretation Eq. (4.12) as a mass effect 
are of general validity at one loop. One can easily check them in other examples, e.g. in 
an effective theory of photons and neutral massive scalars or pseudoscalars. 

At one loop only two scenarios exist. In a massless theory the test particle acquires 
a mass and nothing else. Its speed is modified accordingly. In a massive theory the test 



energies. This can be interpreted as due to a change in the vacuum energy density as 
done in section 2. At high energies, where the mass of the theory should be irrelevant, it 
acquires a mass; precisely the one it would have acquired in the massless theory. Its speed 
is slowed down accordingly. 

From the results presented in this section, only the one corresponding to the Weinberg- 
Salam theory can be straight-forwardly used for phenomenological estimations. The cor- 
rections to the speed of neutrinos are even smaller than the ones of photons reviewed in 
the introduction. 

5. Two-loop QED 

We now turn to the study of the modification of the speed of photons in QED at two 
loops. The calculation has already been performed in the literature for low-energy photons 
[4] [6] . In such a case, the result can be obtained from the Eulcr-Hcisenberg energy density, 
which encodes the fermion box diagram for constant external fields. At high energies, there 
is no possibility to circumvent the use of the full machinery of quantum field theory at finite 
temperature. Therefore, we proceed now to perform that first principles computation. 

For QED, because in the limit in which only massless particles thermalize ^\]^t — 0' 
one can use at two loops all the one-loop formulae and techniques. Thus from Eq. (3.18) 
and Eq. (4.2) one obtains [4] 



particle remains massless at low energies and its speed is slowed down 




f d k 1 

f{q,T,me) = -ie'^n / — — ^ 5{k'^)A{q,k,me) 




(5.1) 



where the integrand, which corresponds to the electron box loop, is given by 




(5.2) 



D{0)D{q-k) D{0)D{q)D{k) D{0)D{k)D{q + k) D{0)D{0)D{k) 
64 mg 64 mg 16 ml{5 q ■ k + 2 ml) 



D{k)D{q)D{k + q) D{0)D{0)D{k)D{q) D{0)D{k)D{q)D{k + q) 



where n is the dimension of space-time and 



D(k) = (p- k)'^ -ml + ir] 



(5.3) 
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The integral in Eq. (5.2) is ultraviolet (UV) convergent. Using Feynman parametrization, 
it is easy to carry out the fermionic momentum integral as well as all the parametric 
ones but the last. This last parametric integral would lead to complicate dilogarithmic 
expressions and, therefore, we prefer to give its result in terms of series. 



A{q, k, me) = ^Yl f(27T3) ^ "^^^^ ^ ^ " 



ml 



{{-ly + 1) (16Z^ + 73/3 ^ ;^Qg;2 ^ ^ 24)] 

(5.4) 

where F is the gamma function. It is clear that the series readily produces, after integration 
upon the thermalized momenta, the low-energy expansion 



2x e^T^ ^ r2(2n)r(2n + l)C(2n + 2) Ari^ + 12?i2 + 5n + 1 f AqT 
« ml) = --^ g r(4n+ l) ^^TT^ 

^ V4050 mj ^ 212625 '^'"J 



(5.5) 

where C stands for the Riemann zeta-function. This gives the expression Eq. (1.3) for the 
leading modification of the speed of low-energy photons. Further higher order terms can 
be obtained without trouble. For arbitrary q the computation is much more involved since 
we need to sum the above series [10]. This is sketched in the Appendix, the basic idea 
being to go back to (5.4) and subdivide it into simpler series which can be summed using 
logs, dilogs and trilogs. The result leads to a high-energy expansion whose first term is 

f{qT » ml) = -^e^T^ In' ^ + O (t' In ^) , (5.6) 
487r^ m^ \ m^J 

leading to 

D \q J rrig 



viqT»ml) = l--{-) In^^ , T<q. (5.7) 



In the process of obtaining this limit, a non-trivial cancellation of In^ takes place. 

The summary of the above computation is that low-energy corrections are of the 
same type as already obtained at one loop order (cf. Eq. (4.9)), but this is not so at 
high energies. Indeed Eq. (5.7) has an IR-sensitive In^ factor which is not present in Eq. 
(4.10). The electron mass plays the role of an IR regulator and the limit me cannot 
be taken. This is consistent with the fact that the m^ — > limit is incompatible with the 
approximations used (which led to H^j^^ = 0). The massless theory cannot be obtained 
from our results and has to be worked out separately. We will do so in the next section 
for a simpler theory. Notice that Eq. (5.6) does not satisfy Eq. (3.13) so that one cannot 
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interpret it as a genuine modification of the speed of fight but rather to one corresponding 
to an induced high-energy mass which runs according to 

X\qT»ml) = ^THn'^. (5.8) 

The high energy correction Eq. (5.7) is never quantitatively more relevant than the low- 
energy correction, Eq. (1.3), as 

Massive theories do not offer better perspectives for obtaining relevant modifications of 
the speed of light at high energies than at low energies. Furthermore only at low energies 
is the modification genuine, i.e. unrelated to an induced mass. 



(5.9) 



6. Two— loop massless gcp'^ 

The results we have obtained so far stress the relevant role played by masses in thermal 
corrections. For QED, me has either compensated for the dimensions carried by the 
temperature or played the role of an IR cutoff. It is, therefore, natural to ask the question 
of what would happen in a theory without any extrinsic mass. This turns out to be a 
subtle problem at two loops, that we will study in the laboratory of massless which 
is substantially simpler than massless QED. 

Let us first recall the one-loop results, coming from the tadpole diagram 

n(i) = , 



24 



where we have used dimensional regularization, which preserves at zero temperature the 
masslessness of the theory. At two loops, Umezawa's two component formalism has to be 
used [18]. This means that Eq. (4.1) is generalized to 



DT{k,k') = 5^^\k - k') 



k'^ + it] (1 _i) \0 -1 



1 




(6.2) 



Fields are accordingly doubled so that the upper component is the physical field and the 
lower, the auxiliary (and always virtual) field which only couples to itself with a quartic 
coupling of opposite sign. Therefore, physical fields only connect to auxiliary ones through 
the off-diagonal thermal propagator. The zero temperature two-loop contributions are 

n(2'*) = , (6.3) 
13 



for the double tadpole, and (n = 4 + 2e) 
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24(167r2)2 



~ ~ — + ci 1 (go - 9 ) 



(6.4) 



for the setting sun diagram. The constant ci depends on the precise definition of the 
regularization scale v. Eq. (6.4) of course vanishes on the mass shell. Since Eq. (6.1) is 
constant, the two-loop generalization of Eq. (3.2) is 



24 



(6.5) 



where we have made explicit the bare character of the coupling constant for the one-loop 
contribution. 

Consider now the diflFerent thermal contributions to the two-loop tadpole diagram. 
With the help of Eq. (6.2), the two-point loop gives a contribution proportional to 




e T — 1 

2 



e T — 1 



so that, recalling 



k'^ + irjj _i \k'^ + ir] 



5ik^)P^ = 



+ m5{k'^) 



(6.6) 



(6.7) 



only the zero temperature contribution is left. Thus the whole thermal contribution comes 
from the other loop. Note also that the Feynman integral 

1 



(27r)" [k'^ + ir^y 



(6.8) 



vanishes in dimensional regularization due to a cancellation of UV and IR divergences. 
As temperature does not affect renormalizability we know that Eq. (6.5) is UV-finite. 
It might, however, be IR-singular on the grounds of being a massless theory which is, 
furthermore, considered on the mass shell. It is, therefore, important to tell IR from UV 
divergences. This we do by introducing where needed the one-loop thermal mass as an IR 
regulator in the two-loop contributions. In particular, Eq. (6.8) becomes 



/ 



(Pk 



{Airy 



In — 5- + C2 



(6.9) 



Here C2 is a constant which depends on the precise value of the IR regulator, which is 
actually arbitrary. Putting everything together, one finds 



n 



(2,t) 



2^2 



1 



1 



a2_ 

48 1 67r2 V e 



+ ln 



- C2 



(6.10) 
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It is worth mentioning here that both in Eq. (6.1) and Eq. (6.10) the thermal loop integral 
has been performed in n = 4 dimensions. It will immediately become clear why there is 
no need to dimensionally regularize it. 

The thermalized setting sun diagram has two real contributions. One corresponding 
to one thermalized particle. It is UV-divergent and IR-finite: 



(6.11) 



Again, the thermal loop integral has been performed in n = 4 dimensions. The other 
contribution, corresponding to two thermalized particles, is both UV and IR-finite. It is 
proportional to the integral 



Jo e -i Jo ey-i J_i J_i 



(6.12) 



^ {r]x + r]'x + fix) + rjrj'xy + yp + - 2^xfj, , 



which reduces to 



Q JO 
ln(l - e" 



dx 



1 



T 

X 1 



+ x In 



+ X 



+ I — - x 1 In 



— X 



2 - 1 



(6.13) 



The high-energy and low-energy limits of I{q,T) are 



(6.14) 



and 



I{q < T) ~ Ag^T 



. n=l 



2 7 3 



where 7 is Euler's constant. Putting together Eq. (6.10), (6.11) and 



(6.15) 



(6.16) 



we get, from Eq. (6.5), 



goT' , g'T^ 1 



E^ = q^ + 



24 



+ 



48 1 67r2 V e 



3 , qT^ ^, qT 
+ ln^ +21n-i^ + c + 



16(27r)^ 



/(g,T). (6.17) 
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Recalling the coupling constant renormalization 

where /i is the renormalization point, we find 

= ,^ + ^ + f in ^ + 21n f ) + ^n,, T) . (6.19) 



24 7687r2 V //^ /^V 25677^ 

Notice that, due to the high-energy and low-energy behavior of /(g,T), Eq. (6.14) and 
(6.15) respectively, its contribution to Eq. (6.19) can be neglected, as it is subdominant. It 
is now also clear why in Eq. (6.5), Eq. (6.10) and Eq. (6.11) there was no need to dimen- 
sionally regularize the thermal Feynman integral: as eventually all UV divergences cancel 
in Eq. (6.19), the 0(e) terms, which one can see that they factorize, do not contribute. 
The /i-dependence in Eq. (6.19) is of course fictitious, so that we can immediately renor- 
malization group improve it by substituting [gT^q^)^/^ for jj?, so that finally (neglecting 
%,r)), 

= (^^ — — ^ . (6.20) 

Thus, the whole two-loop contribution is contained in the running of the coupling constant 
which appears at one loop, 

with A the intrinsic mass scale of the theory. Eq. (6.20) with Eq. (6.21) is the main 
result of this section. The whole two-loop temperature effect is to make the one-loop 
thermally-induced mass run. For very high values of q, 

327r ^ 

{g r^g2)i/3 < A^e^^ , (6.22) 

one starts being sensitive to the Landau pole and perturbation theory breaks down. Below 
that, 

{g T4q2)i/3 ^ ^2 ^ (g_23) 

the running is too little conspicuous to the noticed: the whole temperature effect is a 
constant mass. At much lower values, 

3271- 

(g T^q'^y/^ < A'^e~^^ , (6.24) 
the running coupling vanishes logarithmically and 
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where we have made exphcit that still the conditions corresponding to Eq. (3.5) hold, i.e. 
the kinetic energy still dominates over the mass. From Eq. (3.5), one finds 



v=l + 



< 1 , T « g , 



(6.26) 



which corresponds to the modification of the speed of a massless particle which has acquired 
a logarithmically vanishing mass. As q becomes even lower, q <^ T, an IR divergence 
appears in Eq. (6.19). It requires an IR regulator, which comes in both recalling that 
propagators can be regulated with the one loop thermal mass or recalling that due to the 
induced thermal mass the situation of the pole is shifted. In any case, q^ friezes out at 
gT^ and we find a constant low energy mass 



with the corresponding non-relativistic velocity, Eq. (3.7). 

It is worth stressing again that Eqs. (6.19) or (6.20) are IR sensitive results, as also 
happened, somewhat differently, for the two loop QED results. For gcj)'^ the two-loop 
temperature effects are basically absorbed by the running of the one-loop thermal mass, 
which does so according to the zero temperature one-loop running of the coupling constant. 
There is, thus, no genuine modification of the speed of the massless particle. We surmise 
that this result is of some generality, at least for theories like massless QED which have the 
same asymptotic behavior. They might hold, mutatis mutandis, for asymptotically free 
theories like pure QCD, although more serious IR problems should be expected for them. 

7. Casimir vs. temperature 

We now turn to the study of the propagation of photons in a Casimir vacuum, defined 
by the presence of two infinite, ideal conducting plates. The modifications of observable 
quantities are ultimately related to the careful implementation of boundary conditions. 
In a way, the same is true for thermal effects. There, periodic boundary conditions are 
imposed on imaginary time. We shall keep this analogy present all along our presentation. 

Let us start by comparing Eq. (1.3) ( propagation in a heat bath) with Eq. (1.4) 
( propagation perpendicular to the Casimir plates, cos^ = 1). Clearly both expressions 
are identical with the substitution 2T a~^, but for a sign. This is not surprising as 
temperature corresponds in the imaginary time formalism to periodic boundary conditions 
and Casimir plates to half-periodic ones in one space direction. As said above, one expects 
a correspondence between the formalisms which makes this identity apparent. This was 
already pointed out by Barton [5] . In the same spirit but within a different formalism we 



(6.27) 



9 In 



g2/3T2 
A2 
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want to obtain the analog to Eq. (4.2) for the Casimir vacuum. This will also allow us to 
gain more understanding of what a modified vacuum is. 

The photon propagator for the Casimir boundary conditions Ei = E2 = = on 
= ±1 has been given in ref. [20]. It has the form 

D^''{k,k') = 5^'^\k-k')D^''{k)+d^^\k-k')D^''{k-k3,k's) , (7.1) 

where D^^{k) is the standard free propagator and D^'^{k; ks, fcg) is the gauge-independent 
modification due to the boundaries. The actual form for this second piece of the propagator 
was gotten in ref. [20] in position space 

L»r(^ - y; X3, yz) = - I ^e-'''<^-y^ (-^^ + — - — ^ 

a y 3,</3; ^27r)3 \^ ^ T^{k)) T{k)s\n{aT{k)) (7.2) 

^g-ir(fc)|x3+y3| _ ^iT{k)\a-xz+y3\ _ ^iT(k)\a+xz-y3\ _j_ gjr(fc) Iss+ys | ^ 

where T{k) — k'^ + irj. Since k = {ko,ki,k2, 0), it is clear that momentum is not con- 
served in the 3-direction, as corresponds to the breaking of translational invariance in 
this direction. On the other hand, no energy conservation violation takes place at finite 
temperatures, because the boundary conditions are imposed in imaginary, not real, time. 

Note that Eq. (7.1) includes all the effects due to reflection at the boundaries. As 
mentioned for finite temperatures after Eq. (4.3), one includes in the modified vacuum only 
those effects on the propagation which cannot meaningfully be separated from the standard 
vacuum propagation, i.e. those contributions which also conserve energy and momentum. 
For a Casimir vacuum this means studying the propagation far from the plates. This 
implies to consider a wave packet localized in a region far from the boundaries so that the 
limit fca — > 00 should be taken. From the above expression for the propagator, this limit 
can be taken, obtaining 

D^^ik, ks, k's) - S{k^) \s{k3 + k's) - e'^^'^^^Siks - ki^)] , ka ^ 00 . (7.3) 

In the above limit only two directions are coherently enhanced, one of them forward and 
which thus conserves momentum. This is the contribution which one cannot separate from 
D^^'^ and is, then, included into the Casimir vacuum structure. This leads to the operative 
propagator 



D>^''{k,k') = 5^^\k-k') 



F + ie Y r^(~k) J e-2ir(fc)a _ i 



(7.4) 



where D^^{k) has been written in the Feynman gauge and g^^ — g^'^ + S^^^S'^^. 

The similarity of Eq. (7.4) with Eq. (4.1), with the substitution 2T — > a~^, catches 
the eye, except for the fact that in Eq. (7.4) the exponential has an imaginary exponent 
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and that the tensor in the a-dependent part of Eq. (7.4) is a projector over physical 
hehcities only in the 3-direction. This explains the cos^ dependence of Eq. (1.4), while 
the imaginary exponent implies a Wick rotation when the loop integral is performed in the 
Casimir case which did not happen for the temperature vacuum, leading to an analytic 
continuation (f' — > —(f'. Our results Eq. (5.5) and Eq. (5.6) would, then, become 

« "«) = 4515 (^fei 



and 



leading to 



and 



f{-»ml) = -r^,e^THn- (7.6) 



a 487r^ o 



e 



vi- < mi) = H ( ) > 1 (7.7) 

8100 ^amj ^ ' 



^(^»m^) = l-^(— )'ln^^ <1, (7.8) 

always for propagation perpendicular to the plates. 

The low-energy expression, with v > 1, has been interpreted in section 2 as due to a 
decrease in the vacuum energy density. The high-energy expression, with v < 1, has been 
interpreted in section 5 as due to a high-energy induced mass which, of course, always 
decreases the speed of light. Both results, Eq. (7.7) and Eq. (7.8) required, as discussed 
before, a g ^ 1 as well as a me ^ 1, the condition equivalent to coldness which, here, 
means that the electron Compton wave length should be smaller than the distance between 
plates. These two restrictions not only do not contradict the inequalities that hold for Eq. 
(7.7) and Eq. (7.8) but, rather, ensure that the corrections are small so that perturbation 
theory holds. 

On the contrary, for a cold heat bath (T <S me), the other condition, T -C g, does not 
follow. This, again, shows that a boundary condition in imaginary time does not impose 
kinematical constraints as does a boundary condition in real space. In fact, q <^ T was 
considered in Eq. (6.21). 

Physically, v > 1 in Eq. (7.7) can be interpreted, in harmony with the analysis of 
section 2, in terms of a decrease in the number of modes, recalling that real boundary 
conditions suppress modes. Instead, the presence of temperature adds real particles which 
allow for coherent scattering and, thus, gives v < 1. On the contrary, the result v < 1 in 
Eq. (7.8) can be understood by recalling that in the standard vacuum the contributions 
of the different modes balance in such a way that photons remain massless but that now, 
by suppressing some modes, this balance is broken and a high energy mass appears. This 
certainly slows down the propagation of photons. 

It is worth mentioning that examples of repulsive Casimir effects are known. However, 
a detail analysis [21] show that the energy density remains lower and, consequently, the 
speed of light is faster than c. 
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For gravity, the result > 1 [2] [8] can be interpreted similarly as a decrease in the 
number of modes, and v < 1 as the appearance of a mass, following the same line of 
thought as done above for the Casimir vacuum. 



8. Conclusions 



We have studied, up to two loops, the modifications of the speed of massless particles 
in a cold heat bath and in other modified vacua, for massive and massless theories, for 
high and low energies. We have found a unifying expression for the low energy change in 
the averaged speed of a light valid for all vacua: 

where p is the energy density and where, if the vacuum is a gravitational one, one a has to 
be substituted by rn^GN. It follows automatically that if the vacuum has a lower energy 
density than the standard vacuum, p < and v > 1, and viceversa. 

An enigmatic consequence of assuming full generality of the equation we have just 
discussed takes place in cosmology. There, open, closed and critical universes are distin- 
guished by their total energy content. It follows from Eq. (8.1) that the critical universe 
is singularized as low-energy photons travel in it with speed c, due to the cancellation of 
all vacuum energy contributions. 

At one loop the whole temperature effect in the case of a massless theory is just to 
give a mass. For a massive theory it leads to a genuine change in the speed of the particle 
at low energies, but suppressed by (^) . It parallels the two-loop results in Eq. (8.1). 
At high energies, it leads to a change in the speed of the particle as caused by a thermal 
mass, which only plays a role at high energies. The effect is suppressed by (|^) . 

At two loops the whole temperature effect in the case of the massless theory chosen, 
gcp"^, is to let the thermally induced one-loop mass run, through the running of the coupling 
constant. It is an IR-sensitive result which requires an IR cutoff, which we have taken 
to be the one loop-thermal mass. All the changes in the speed of the massless particle 
are due to the now running thermal mass. For the massive theory chosen, QED, the two- 
loop low-energy result reproduces the known modifications, and is suppressed by the same 
factor as at one loop. The high-energy behavior, however, differs from the one obtained 
at one loop by a squared logarithmic correction which, furthermore, depends on the mass 
of the theory in an IR sensitive way. Nevertheless, it can be interpreted, like at one loop, 
as due to a now running thermal mass. 

Finally, the modification of the speed of light due to Casimir plates has been compared 
and related to the one due to temperature. A field theoretical explanation in terms of 
modes suggests the following physical picture of why photons move faster between plates 
than in a normal vacuum, in contrast to what happens in a heated vacuum. Modifications 
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of the vacuum that populate it with real or virtual particles introduce coherent scattering 
which reduces the speed of massless particles. The opposite situation, e.g. Casimir vacuum, 
corresponds to modifications of the vacuum that eliminate virtual modes and, consequently, 
their would-be scattering. The speed of massless particles is, then, increased. 

The most concise summary of our results states that only low-energy massless par- 
ticles corresponding to a massive theory show genuine modifications of their speed while 
remaining massless. All other modifications are mass-related, or running mass-related. 
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Appendix 



In section 5, we needed to calculate 



j{q,i,me) = -■ 

where 



m ^ I f^f^V ((-1)' ^ 1) (i + (ITWTi) ) ■ 

First, we define 

R{y) = S{y) + S{-y) , (8.4) 
with S containing all the powers in the series 

^M-Efg^^^"V(l^ (,^g,^,) ). (8.5) 
This sum can be splitted into three pieces, 

S{y) = \ + bS2{y) + 2SM) ■ (8-6) 

y 

The first two are easily carried out in terms of logarithms. The third requires changes of 
variables to get it ready for resumation in terms of trilogarithms. These partial results are 

2y3 



= -2/ - - ^ - 2i^y{l-y)\n{i^ + + (1 - 2y) ln\iy/^ + , 

(8.7) 

c? - r(n-i)r(n-2) r 

^l + ly^- Iv"" + (1 + 2y)Wy{l-y)ln{i^ + y/T^) + (-^ + 2y) ln\i^ 

+ aA^) , 

(8.8) 

= I - f - Y + - 18y)ln(^Vy + yr^) 

+ ( - ^ - 2y - 2y2)ln2(z^+ y^) + l^ln3(i^ + ^13^) (8.9) 

+ Ay In'ii^ + ^/T^) H-2iy/l-y^ + 2y) 

- 4y ln(iVy + Vl - y)L2 [iV^-V + ^Vvf] 

+ 2yL3 [(aA^ + iVy] - 2yC(3) , 
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where we have used the notation 



Lniv) = E F • (8-10) 
fe=l 



Collecting this partial results we obtain 



II- --2 14 ,3 



+ (4 - 8y)V2/(l-y)ln(iVy + 

- 8yln2(i^+ ^Jl^) \n{-2i.Jy{l-y) + 2y) 

- 8y ln(z^ + v'r^)L2 [(iVy + 



It is then easy to compute the limits used in the text, 



R{y) ~ -21n"y 

j/->oo 

88 2 416 4 

~ — V -\ y + . . . 

y^Q 135^ 3375^ 



(8.12) 
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